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Electronic transport properties of Au nano-structure are 
investigated using both experimental and theoretical analysis. 
Experimentally, stable Au nanowires were created using me- 
chanically controllable break junction in air, and simultane- 
ous current-voltage (I-V) and differential conductance 51 /SV 
data were measured. The atomic device scale structures are 
mechanically very stable up to bias voltage 14 ~ 0.6V and 
have a life time of a few minutes. Facilitated by a shape 
function data analysis technique which finger-prints electronic 
properties of the atomic device, our data show clearly differ- 
ential conductance fluctuations with an amplitude > 1% at 
room temperature, and a nonlinear I-V characteristics. To 
understand the transport features of these atomic scale con- 
ductors, we carried out ab initio calculations on various Au 
atomic wires. The theoretical results demonstrate that trans- 
port properties of these systems crucially depend on the elec- 
tronic properties of the scattering region, the leads, and most 
importantly the interaction of the scattering region with the 
leads. For ideal, clean Au contacts, the theoretical results in- 
dicate a linear I-V behavior for bias voltage 14 < 0.5V. When 
sulfur impurities exist at the contact junction, nonlinear I-V 
curves emerge due to a tunnelling barrier established in the 
presence of the S atom. The most striking observation is that 
even a single S atom can cause a qualitative change of the 
I-V curve from linear to nonlinear. A quantitatively favor- 
able comparison between experimental data and theoretical 
results is obtained. We also report other results concerning 
quantum transport through Au atomic contacts. 

I. INTRODUCTION 

Electron transport through atomic nano-contacts has 
been an active research area for a decade both experi- 
mentally and theoretically. The scientific interest of these 
systems is largely driven by their peculiar electronic and 
transport behavior. The atomic nano-contacts are struc- 
tures with low atomic coordination number and, as a 
result, can behave very differently from the bulk coun- 
terpart. From a practical point of view, understand- 
ing the novel electronic and structural properties of the 
atomic nano-contacts is an important step towards nano- 
device fabrication and characterization. The first set of 
experiments on nano-contacts focused on their zero bias 
conducta^e (G) using Scanning Tunnelling Microscopy 
(STM)yij^echanically Controllable Junction 
(MCBJ)Era, and Relay Contacts (RC)BE3. In these 
experiments a narrow constriction with a few atoms at 
the cross section is formed. As the electrodes are pulled 
apart, G is measured and found to change discontinu- 
ously forming plateaus with values close to n x Gq; where 
n is an integer and Go = 2e^//i ~ Ikh conduc- 



tance quanta. Pioneering experimentsollj clearly showed 
the correlation between conductance jumps and mechan- 
ical properties in the nano.-contacts. These results con- 
firmed earlier predictionscJ that the conductance varia- 
tions are due to abrupt changes of nano-structure cross 
section as a function of wire elongation. Extensive theo- 
retical investigations on nano-structures have been pub- 
lished recently to analyze these systems. One major 
theory focus is to calculate the zero bias conductance 
through a ballistic quantum point contact j—, These cal- 
culations start by assuming various contactES^E^I geome- 
try, or by using more realistic atomic positions derived 
from molecular dynamics simulations. The potential of 
the constriCrticp, and/or interaction Hamiltonian is then 
constructedcTH from which the zero bias transport co- 
efficients are evaluated. While different levels of approxi- 
mations were employed in these theoretical analysis, den- 
sity functionaJ—thepty based ah initio analysis have also 
been reportedc3^E£lE3 which provide self-consistent calcu- 
lations of atomic nano-contacts. 

While zero bias transport coefficients have received a 
great deal of attention, one must go beyond this limit 
to understand the full nonlinear current-voltage (I-V) 
characteristics of the nano-contacts, as this information 
is essential for the understanding of real device opera- 
tion. For example, due to the small cross section of these 
systems, they are exposed to substantial current density 
~ IQ^A/cm?' which may result in atomic rearrangement. 
It is also expected that electron-electron (e-e) interaction 
is enhanced due to the strong lateral confinement pos- 
sibly leading to Luttinger liquid behavior for the quasi 
one dimensional atomic wires. Separating these differ- 
ent effects is experimentally challenging due to the many 
variables which can affect the results. This is probably 
the origin of the existing controversy in explaining the ex- 
perimental otjserved non-linear I-V curves of the atomic 
scale wiresE^T^. From a theoretical point of view, this is 
also a challenging problem: so far only two computation- 
ally accurate techniques exist which can treat systems 
with npep-boundaries out of equilibrium ducj-tn-external 
biafS'L3'L3. In the approach of Taylor et. alE3'cB, realis- 
tic atomic leads can be treated and the problem is solved 
self consistently within LDA. Therefore, the leads, the 
device (scattering region) and their couplings are incor- 
porated without any preconditioned parameters. 

To further shed light on the physics of quantum trans- 
port at molecular scale, we report in this paper our in- 
vestigation on transport properties of Au nanostructures 
both experimentally and theoretically. Experimentally, 
we created stable Au nanowires using the mechanically 
controllable break junction technique in air, and we si- 
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multaneously measure the I-V curve and the differential 
conductance 51 /5V . We found that our atomic scale Au 
nano-contacts are mechanically very stable up to bias 
voltage Vb ~ Q.&V and have a life time of a few minutes 
which is adequate for our measurements. As we are in- 
terested in features due to electronic degrees of freedom 
of the nano-contacts, careful data analysis is needed be- 
cause transport data can be affected by many factors. 
By defining the shape function, S = -^j-^ which finger- 
prints electronic properties of the atomic device, our data 
clearly shows differential conductance fluctuations with 
an amplitude > 1% at room temperature, and a non- 
linear I-V characteristics. To understand these trans- 
port features, we carried out ab initio calculations on 
various Au atomic wires bonded with atomic Au elec- 
trodes using the first principles technique of Ref. |7 58 
Our calculations show that pure and perfect Au nano- 
contacts do not give the nonlinear I-V curves as mea- 
sured in the experiments. However when Sulfur impu- 
rities are present near the wire-electrode contact region, 
the nonlinear I-V curves emerge due to the tunnelling 
barrier provided by the impurity atoms. The most strik- 
ing observation is that even a single S atom can cause a 
qualitative change of the TV curve from linear to non- 
linear. A quantitatively favorable comparison between 
experimental data and theory results is then obtained. 
Our combined theory-experiment investigation allows us 
to conclude that transport through Au atomic wires is 
strongly affected by the properties of the wire-electrode 
contacts. 

The rest of the paper is organized as follows. In the 
next section, experimental measurement and results are 
presented. Section III presents the theoretical results 
while Section IV discusses the transmission coefficients 
in more detail. We also discuss and compare previous 
works with ours in section V, followed by a conclusion. 



II. EXPERIMENTAL RESULTS 

Atomic scale gold contacts were formed with a mechan- 
ically controllable break junction in air at room tempera- 
ture. Once suitably stable atomic junctions were formed, 
a slowly varying bias voltage was applied (typically a 
O.lHz triangle wave, 2Vp-p) along with a small mod- 
ulation voltage (typically AOkHz, AmVrms) across the 
contact and a load resistor of 3Kil. Current (/) and 
differential conductance (Sl/dV) were measured with an 
TV preamplifier and a lock-in amplifier. The experimen- 
tal set up is shown in Figjl|-a. A typical measurement 
through a stable Au nano-contact is presented in Fig.^ 
b. We show our data for the differential conductance and 
the DC conductance (G = I/V). In the inset of Fig|l]-b, 
we also display a typical I-V measurement. This data was 
taken over the course of a 5 seconds voltage sweep from 
positive (dark lines) to negative (grey lines) bias voltage 
(V). Both polarities share a common overall shape, but 



seem to vary significantly in the details of their conduc- 
tance behavior. However, for each polarity, one notices 
that there seem to be similar details present in the DC 
and differential conductance. 

An important issue in order to understand these results 
is the separation of effects due to atomic rearrangement 
in the nano structure from electronic properties. We ad- 
dress this problem by writing, very generally: 



IiX,V)=giX)f{X,V) 



(1) 



where the variable X symbolizes the effects of atomic 
structure, but no explicit knowledge of X is required. We 
define g{X) as the unbiased conductance and f{X, V) is 
the normalized functional form of the voltage dependence 
of the current. Therefore, in the zero bias limit V ^ 0, 
f{X, V) = V. We further define a new quantity called 
the "shape function" {S), 
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Fig.^a plots S corresponding to the data of FigJ^b, for 
both positive and negative bias voltages. The similarity 
between the curves of the shape functions for both polar- 
ities, including fine details, is in striking contrast to the 
easily distinguishable conductance plots of Figj^b. By 
its definition (Eq. ||), S depends only on /, the functional 
form of I{X,V), and is independent from the zero- bias 
conductance g{X). This fact and the usefulness of the 
shape function S can be seen by considering the follow- 
ing "Gedankcn experiment". Let's assume that we mea- 
sured / and dl /SV through a variable, ohmic potentiome- 
ter as a function of the applied voltage. Suppose there 
were drastic changes in the temperature during the mea- 
surement, and some troublemaker stochastically turned 
the knob of the variable potentiometer without telling 
anyone. Glancing at the measurements of I{X, V) and 
51 /6V alone, one might wrongfully conclude that the po- 
tentiometer was exhibiting a nonlinear behavior. How- 
ever, a plot of S would show that S{V) — 1 (from Eq.^), 
which would allow us to deduce that the I-V curve was 
actually linear and thus in fact ohmic. One would also 
conclude that the origin of the apparently nonlinear be- 
havior was due to a change in the unbiased conductance 
g{X), rather than due to a true nonlinearity in voltage. 

Returning to the experimental results shown in Fig.^ 
a, we note the overlap of the two curves of S for positive 
and negative bias voltages; this points to a similar func- 
tional form for both bias polarities. Because of this, we 
deduce that for this case f{X, V), at most, depends very 
weakly on X (especially at bias voltage Vb < OATV). We 

can therefore write f{X, V) ~ foiV) e-f , and the 
unbiased conductance g{X) — ^j^y)"* ■ Hence, as an ex- 
perimental voltage sweep typically takes 5 seconds, time 
dependent atomic rearrangements (changes in X) man- 
ifest themselves as stochastic variations of g{X). Other 
details of this data analysis technique and further dis- 
cussion on the shape function are not within the scope of 
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this paper and can be found elsewhereEj. 

The normahzed functional form of the voltage depen- 
dence of fo{V) for the measured data is shown in Fig.||- 
b. We note that the curves for both polarities appear 
on the top of each other and they are indistinguish- 
able. However, as expected, the "trouble maker" in our 
"Gedanken" experiment shows up as fluctuations in g{X) 
shown in Fig.^c, indicated by the fluctuations and by the 
fact that positive and negative bias give different traces 
of g{X). Note that during the course of our bias sweep 
from initial bias voltage to some bias value V = Vi, the 
atomic structure has been fluctuating and changed from 
what we started with to something unknown. But if we 
could freeze the structure at that moment and re-measure 
the I — V curve for the fixed structure, g{X) would be 
its zero bias conductance. In other words. Fig. |^c shows 
a parametric plot of what the conductance g{X) would 
be at zero bias at the point in time when this voltage 
was measured experimentally; hence g{X) could also be 
viewed as a function of time in this figure. Our analysis 
show that the fluctuations of g{X) which are less than 
5% peak-peak, are attributed to changes in the atomic 
configuration of the junction. We have also found that 
g{X) exhibits no time correlation and its Fourier Trans- 
form has a 1/f frequency behavior. 

In Fig.^d we plot the normalized differential conduc- 
tance, Sf /SV. Examining the curves in Fig.^d, we can 
see the same subtle fluctuations in Sf / SV as those found 
in S (Fig||-a). Both plots show broad and fine details 
that are uncorrelated to the fluctuations of the unbiased 
conductance g{X) (Fig|2|-c). These plots reveal that it is 
the fluctuations of the g{X) which dominate the fluctu- 
ation features in the measured I{X, V) and SI/6V. 

Over different voltage ranges (0.2 — 0.35V, 04, — 
0.35V, 0.1 - 0.5V), Sf/SV has a high correlationEJ of 
0.99 between both polarities. Over these same ranges, 
g{X) had a weak correlation (—0.5, 0.25, —0.2, respec- 
tively) between both polarities which changed drastically 
depending on the selected voltage range. We have also 
calculated the correlations by shifting the voltage of one 
polarity with respect to the other by ±10mV (the scale 
of the flne details) in increments of ImV. Over all three 
ranges, the correlation for Sf/SV had a local maximum 
for OV shift. No correlation extrema were found in the 
case of g{X). This quantifies how similar the details in 
both polarities of 6f /5V are to each other, in contrast to 
the more easily distinguishable curves for g{X). We thus 
conclude that the "wiggles" observed in ^ are electronic 
in nature and they are not variations in g{X). 

The "wiggles" (magnified in the inset of figure |^d) 
of Sf/SV may actually be much more pronounced than 
indicated in these plots. They are smeared out by un- 
avoidable experimental constraints. Since we must add 
a modulation signal to make our lock-in measurement of 
51 /SV, we end up averaging over a range of V. This leads 
to broadening and decrease in amplitude of these fine 



detaila^j'Eil. This unavoidable averaging artifact keeps us 
from making a more precise measurement of the shape, 
amplitude and voltage characteristics of this fine struc- 
ture, but the true features should be sharper and more 
pronounced than they are measured to be. From our 
data, we observe "wiggles" on the voltage scale < lOmV 
and with amplitude ~ 1% of the signal. The width of 
these fluctuations is comparable to the modulation volt- 
age, AmVrms ^ IQinVp-p] hence, it is possible to have 
features on a voltage scale < lOmV with amplitude which 
is orders of magnitude larger than these 1% fluctuations. 

If one were to merely examine our /(V) measurements, 
similar results have been , observed in air and at low 
temperatures with RCOGj and STM configurationsE3E3. 
The fine details exposed by the analysis presented above 
have not been discussed in the literature as they tend 
to be hidden by the fluctuations in the unbiased con- 
ductance g{X) (see Fig.||-c). The data presented here 
represents the behavior of one junction. However, we em- 
phasize that the general features presented in Fig.|| are 
experimentally found to be independent from the zero 
bias conductance g{V = 0) for values between 1-10 Go- 
The fine details of Sf/SV are reproducible for a given 
stable junction, but the speciflc details change for differ- 
ent junctions. Although there seems to be a voltage scale 
associated with these details (< lOmV), Fourier analy- 
sis does not indicate any strong periodicity. It should 
be noted that we have imposed a selection rule on the 
junction type by studying device conflgurations that are 
stable on the time scale of minutes. We have also ob- 
served junctiona^Lhat exhibit linear behavior {S = 1) 
as have othertSoO. These junctions however are not 
stable over this long time scale. In comparison to their 
nonlinear counterparts, linear junctions have smaller fluc- 
tuations in g{X). 

In the following we will discuss the physical origin of 
the observed /(V) characteristics. In this aspect we will 
investigate the electronic effects rather than the struc- 
tural parameters which will be assumed static. This will 
allow us to gain valuable insight into the voltage depen- 
dent conduction properties of nanoscale electrical con- 
tacts. We will thus compare our modeling to the nor- 
malized functional form of the current, /o(V), shown in 
Fig.^b, rather than the original I-V curve presented in 
the inset of Fig|l]-b. We will also explain the origin of the 
fluctuations in the normalized differential conductance 
shown in Fig.||-d as well as the effects of temperature 
and modulation signal on their amplitude. 

III. AB INITIO ANALYSIS OF THE I-V 
CHARACTERISTICS 

To provide a theoretical understanding of the experi- 
mental data presented above, we have calculated the I-V 
characteristics of Au nano-contacts self consistently by 
combining the density functional theory and the Keldysh 
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nonequilibrium Green's functions. The method is based 
on the newly developed ab initio approach to treat open 
electronic systems under finite bias. For technical de- 
tails of this metlia|d we refer interested readers to the 
original papersEZlH. Very brieflyp-jMXK-analysis uses an 
s, p, d real space LCAO basis setEO'EsS and the atomic 
cores are defined by the. standard nonlocal norm con- 
serving pseudopotentialo The density matrix of the de- 
vice is constructed via Keldysh nonequilibrium Green's 
functions, and the external bias Vh provides the elec- 
trostatic boundary conditions for the Hartree poten- 
tial which is solved in a three dimensional real space 
grid. Once the density matrix is obtained, the Kohn- 
Sham effective potential Ve//(r; VJ,), which includes con- 
tributions from Hartree, exchange, correlation and the 
atomic core, is calculated. This process is iterated un- 
til numerical convergence of the self-consistent density 
matrix is achieved. In this way, we obtain the bias 
dependent self-consistent-|affective potential T4//(r;V6), 
from which we calculateEllGa the transmission coefficient 
T{E, Vb) = T{E, [Veff{r, Vt)]), where E is the scattering 
electron energy and T is a function of bias Vb through its 
fimctional dependence on Ve//(r; Vb). 

In our analysis, the scattering states are defined for 
energy ranges between the left and right chemical poten- 
tials fiL and fiR, respectively. To solve for these states, 
at a given energy—E, we solve an inverse energy band 
structure problemE^I. We then group all states as left 
and right propagating states depending on their group 
velocity. For a scattering state coming from the left lead, 

^f^" should start as a right propagating state and it 
gets refiected back as a left propagating state t/)^'" with 
refiection coefficient r^^'^n in the left lead, and trans- 
mitted into the right lead as a right propagating state 

. . . . 

4>fj^ with transmission coefficient i^'"' "■ In the nu- 
merics, the scattering states are represented as a linear 
combination of atomic orbitals inside the device region. 
This allows us to write, for example, a left scattering 
state as: 

{^^L + ^f^r^"'^" inside left lead 
7/;^" inside device 

ri?/"^ " mside right lead 

A scattering state in the right lead can be written in a 
similar fashion. For a symmetric two probe device, the 
total transmission frona^he left lead is identical to the 
one from the right leadE3. To calculate the total current 
inside the device at a given bias voltage VJ, applied to the 
right lead, we use: 

/(^b) = -^y dET{E,Vb)[fL{E,iJiL= IJio) 

-.fRiE,fiR^^io + eVb)] (3) 

where /^(fl) is the Fermi function on the left (right) 
lead evaluated at temperature T = OK unless otherwise 



stated. In addition to predicting the overall transport 
properties of a device, our formalism enables us to study 
transmission through each incoming Bloch state of the 
leads separately, therefore allowing us to separate effects 
due to the leads and due to the scattering region. We 
have used this formalism to calculate I-V characteristics 
of structurally different Au nano-contacts and compared 
them with the experimental results described above. 

A. Perfect Au nano-contacts 

In a first attempt to model our experiments, we 
calculated the I-V characteristics of four Au atoms 
(= "molecule" ) in contact with Au(lOO) leads. The struc- 
ture of the atomic device is illustrated in Fig.^a. The 
scattering region is bonded by two semi-infinite Au leads 
which extend to electron reservoirs at ±cxd where bias 
voltage is applied and current is collected. The device 
scattering region, indicated by D, is described by three 
Au layers from the left lead, the four Au atoms in a chain, 
and two layers of Au from the right lead. We have also 
increased the two Au layers on the right side of the chain 
to four to ensure that convergence is reached with respect 
to the screening length. We note that long and thin-gpld 
necks have directly been observed experimentallyEEl.In 
this structure, the registry of the atomic chain with re- 
spect to lead surface layer can be different. The most 
common structures, which we analyze in this work are: 
the hollow site, where the atomic chain faces the vacant 
position in the lead layer as shown in the upper left in- 
set of Figj^B; and the top site, where the atomic chain 
and an atom from lead surface layer face each other as 
illustrated in the upper right inset of this figure. 

Our calculations show that in all cases charge transfer 
between the Au chain and the leads is not important, be- 
ing only ^ 0.07—0.1 electrons per atom to the Au chain at 
different bias voltages. This corresponds to less than 1% 
difference in electron population per atom and therefore 
does not play any significant role in the I-V characteris- 
tics of Au contacts. This is in contrast to a binary atomic 
system such as carbon chains between Al(lOO) lcadsE3. In 
addition, solving for the energy eigenvalues of the four 
atom chain gives a HOMO-LUMO gap of 0.68eV, indi- 
cating that the molecule eigenstates should only have a 
secondary effect on the transport properties. Therefore, 
the major effect on the I-V characteristics is due to the 
character of Bloch states in the leads and their couplings 
to the molecule at the chain-lead interface. 

In Figj^A, we show the results for a system with atoms 
in the hollow site. At small bias Vb, we note that current 
is a linear function of Vb with a slope G 0.94Go. The 
linear function suggests that T{E, V) — To with a weak 
voltage dependence. In this regard, our self-consistent 
calculation gives-apcesult apparently similar to previous 
theoretical worknJtj in the low bias regime. However, we 
will show later, by addressing the origin of this "perfect 
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linearity" , that the physical picture of a bias-independent 
transmission coefficient is not valid even for such a sim- 
ple chain. We also note that the linear I-V characteristics 
observed in these systems do not agree with our experi- 
mental data. 

A major feature of FigJ^A is the huge plateau at 
Vb = 0.5V — 0.9V, as well as the fine structures (or 
sometimes negative differential resistance) observed for 
larger voltages. In the lower inset of Fig.^B, we show the 
band structure of the Au(lOO) lead along the z— direction 
(transport direction). Even-though at a given energy E 
there are many Bloch states which are potential candi- 
dates for transporting current, our investigation found 
that for E < 0.5eV there is only one state that is ac- 
tually conducting (presented by a continuous line in the 
inset). Once this state is terminated at E 0.5eV, the 
current is saturated resulting in a large plateau until new 
conducting states emerge at higher bias voltages. For 
E > 0.9V, transport properties are more complex since 
more states contribute to transmission. Under these cir- 
cumstances, band crossing occurs more frequently and 
T{E, Vb) changes over small ranges of bias leading to the 
small structures seen in the I-V characteristics of perfect 
Au contacts. In fact these variations in T{E, Vb) are the 
origin of the fluctuations in the normalized differential 
conductance shown in Fig.||-d. They thus need to be at- 
tributed to the effects of the leads' band structure. In the 
lower inset of Fig.^A, we plot the theoretical SI/SV. The 
magnitude of these conductance fluctuations is of the or- 
der of 60% at zero temperature which is much larger than 
the experimental finding. However, these fluctuations are 
reduced to 15% if the current is calculated using Eq.|| at 
a temperature T = 300K (the temperature of our exper- 
iments). In addition, the fluctuations further decrease 
to ~ 1% when current is averaged over the experimental 
modulation voltage range of AmV, completely consistent 
with our experimental results of Fig.^d. 

The effect of site registry is studied by placing the end- 
atom of the Au chain at the top site of the leads. In 
this situation the chain atoms are facing one atom of the 
leads' surface layer|S-[Xhis analysis is quite important, be- 
cause it was shownEacZl that atoms at the junction change 
registry from hollow to top sites resulting in bundle for- 
mation just before the nano structure breaks. The I-V 
characteristics of these systems (shown in Fig.^B) are 
similar to the previous results, therefore no change is ob- 
served in the transport properties for the top site registry. 

To further investigate the huge and peculiar plateau 
occurring at Vb = 0.5V — 0.9V, we have plotted in Fig.^ 
the charge density at a given energy E, p{E,x,z) = 
J dyp{E,x,y,z). We see clearly that at zero bias Vb — 
OV , the device property turns from a perfect conductor at 
E = OeV (Fig.|-a) to an insulator &i E = O.QSeV (Fig.|- 
c) due to the termination of the conducting state. Here 
we interpret the charge concentration as the effective 
bonding strength, or conductance probability. Applying 
a bias voltage to the system drives it out of equilibrium, 
and at Vb = 0.68^ and E = O.Qd>eV, the charge in the 



molecule redistributes, but the bonding is still very weak 
as shown in Figjij-d, with some molecular regions having 
zero charge and resulting in the large plateau observed 
in our I ~V curve of FigJ^A,B. This effectively demon- 
strates the importance of both the energy and the voltage 
dependence of the transmission coefficient T{E, V). This 
point will be discussed in more detail in section V. A fur- 
ther point to notice is the difference of zero bias conduc- 
tance, G ~ O.SGq for the top-site device and G ~ 0.9400 
for the hollow-site device. This difference is due to a 
change in the coupling between the chain end-atoms and 
the surface of the leads. To ensure the same nearest 
neighbor separation distance for both cases, we end up 
with four nearest neighbors for the hollow site registry 
and only one nearest neighbor for the top site registry. 
Under these circumstances, the hollow site has a better 
coupling to the chain and hence a larger conductance. 

A pure and perfect Au nano-contact, as studied in this 
section, shows rich and interesting transport properties. 
It also gives a good understanding of the origin of the 
observed differential conductance fluctuations as due to 
coupling of the Au-chain to the leads' band structure. 
However, it has a linear I-V curve at small Vb < 0.5V, 
rather than the experimentally observed nonlinear I-V 
characteristics. How can a nano-contact produce a non- 
linear I-V curve such as that of Fig.(|^b) ? The simplest 
possibility to observe such a phenomenon is to have a 
tunnelling barrier at the molecule-lead junction whose 
effect gradually collapses as a function of an increasing 
bias voltage. Indeed, recent experimental findings indi- 
cate that perfectly linear I-V characteristics were repro- 
ducibly -found in gold-gold nano-contacts in ultra high 
vacuumEJ, while nonlinear effects emerge when the ex- 
periment was performed in the air. This suggests that 
impurities play an important factor. Therefore, one has 
to go one step further and investigate the effect of impu- 
rity and disorder at the Au junction. 

B. Au nano-contacts with S impurity 

To simulate the effect of an impurity at the contact, we 
have replaced one of the Au atoms at the interface layer 
with a sulphur atom. This is presented in the inset of 
Fig.^. The choice of sulphur is motivated by the fact that 
in our experimental labs located in down-town Montreal, 
sulphur is a non-negligible airborne pollutant (diesel ex- 
hausts); sulphur atoms bond actively with Au. We also 
note that the band-width of sulphur, ~ lOeU, is much 
higher than that of Au (~ leU), thus a tunnelling bar- 
rier is expected to be provided by the presence of S atoms. 
In this system, charge transfer to the atomic chain is still 
small, and thus inadequate in explaining the experimen- 
tally observed TV characteristics. We note that the S 
atom suffers from an electron deficiency ~ 4%. Also due 
to the presence of the S atom, the coupling of the Bloch 
states in the leads to the device scattering region is quite 
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different as compared to the mono-atomic gold structure. 
When the S atom is present, our calculations found that 
all Bloch states are coupled to the scattering region, with 
the highest transmitting mode still being the one corre- 
sponding to the conducting mode of the mono-atomic 
gold system. 

The TV characteristics of the S doped Au nano- 
contacts is shown in Fig.^. We note that the TV curve 
for voltages up to 0.5V is very similar to the experimental 
values with nonlinearity onset at non zero bias voltage. 
We also note that the huge current plateau of the pure 
Au device has now diminished because more states con- 
tribute to electronic transport. To compare these results 
with the experimental measurement, we have used the 
normalized functional form of the voltage, fo{V) shown 
in Fi^|-b, and multiplied it by the simulated zero bias 
conductance (0.67Go)- The result is shown as open cir- 
cles in Fig.^. The qualitative and quantitative agree- 
ment between the theoretical and experimental results is 
rather encouraging. In fact, this is the first time that 
experimental nonlinear I-V characteristics could be com- 
pared so well and so directly with ab initio self consistent 
calculations. It is also a very surprising result because 
a single S impurity can qualitatively alter transport in 
these nano-contacts from linear to nonlinear. 

The TV curve in Fig.|| still shows the small features 
similar to those found in pure and perfect Au contacts 
(presented in Fig.|^) . These fine details of the calculated 
TV curve would result in differential conductance fluc- 
tuations similar to the ones shown in the experimental 
data of Fig.^d and the pure Au device of Figj^b. How- 
ever, these fine features as calculated are wider and occur 
at higher bias voltages than the experimentally observed 
ones. The former can be attributed to the zero temper- 
ature we used in our calculation. The absence of these 
fine features at smaller voltages is attributed to the small 
size of the leads used in our theoretical modeling: close 
to Ep there are just a few states so that abrupt variation 
of T{E, Vb) at smaller 14 is less probable, resulting in a 
smoother I-V curve at low Vb- 

To understand other possible factors which can affect 
I-V characteristics of nano-contacts we have studied the 
effect of a larger number of S impurities, disorder, and 
their combined effects. The results of these calculations 
are presented in the following subsection. 

C. Contacts with several impurities and disorder 

Including more S impurities at the contact enhances 
the tunnelling barrier and may give rise to a smaller cur- 
rent with a larger nonlinear behavior. The result with 
replacing two Au atoms at the interface by S atoms is 
plotted in Fig.||-A. Indeed, as expected the TV curve for 
this system shows a larger nonlinear character. The non- 
linearity starts at Vb ^ 0.17eV. It is actually more non- 
linear than that of the experimental data; this is mainly 



due to the high concentration of S at the interface. There- 
fore, an experiment which dopes an Au contact with more 
S impurity should show an enhancement of the I-V non- 
linearity. This also gives a possible explanation of why 
experimentally [the nonlinear I-V fitting parameters are 
not universaEj'Efl: they dramatically depend on the con- 
tact structure as well as the impurity concentration. 

Studying the effect of disorder in nano-contacts is an- 
other important problem. To investigate this effect, we 
have randomized the contact layer in the left lead of a 
pure Au nano-contact, as shown in the inset of Fig.^B. 
The distribution of disorder leads to a smaller distance 
between the contact layer and the Au chain, resulting in a 
better coupling. The current through this device is larger 
than that of the ideal contact and is shown in Fig.^B. 
For this device, the slope of the current at zero bias is 
G - 0.915Go, and it slightly increases to G - 0.92Go 
at Vb ~ Q.2V. The I-V curve shows very weak nonlinear 
characteristics. This suggests that disorder alone may 
create a tunnelling barrier which is overcome through the 
application of a bias voltage. However, to observe the ef- 
fect, there need to be conducting states in the scattering 
region. For our pure Au device, only one single state 
is conducting and the maximum zero-bias conductance 
is G = Go- Therefore, at G ^ 0.915Go, the conducting 
channel is already open to near its maximum at zero bias, 
hence it cannot be further enhanced in any significant 
way by applying a bias. The results in Figj6|-B show that 
disorder is an important factor that allows more Bloch 
states in the leads to couple with the scattering region 
and contribute to transport properties. This is clearly 
seen when we notice that the huge current plateau ob- 
served in Fig.^ essentially vanishes in disordered device. 

Combining the effect of disorder and S impurity is also 
crucial. The latter enhances the tunnelling barrier and 
the former may enhance the coupling of the device to the 
leads. We have used the disordered structure studied in 
the last paragraph and replaced one of the Au atoms at 
the contact layer with an S atom, as shown in the insets 
of Fig.^C. The zero bias conductance for this device is 
G ^ 0.833Go. This number is larger than the one with 
only a tunnelling barrier (ideal contact with S impurity), 
but is smaller than the one with only a disordered con- 
tact (which has better coupling). The final conductance 
is due to a competition between both effects. The TV 
curve for this device is presented in Fig.^C, it shows very 
weak nonlinear I-V characteristics: analysis of our data 
show that the slope reaches G ~ 0.85Go at Vb ~ 0.18V. 
Therefore, there is still a weak nonlinear behavior, but it 
is small due to the fast saturation of conducting channel 
in the device. 

From these results we can conclude that the TV charac- 
teristics of an atomic junction is a complex phenomenon 
in which the leads, eigenstates of the scattering region, 
as well as impurity and disorder play major roles. In 
particular, a tunnelling barrier created by impurities can 
result in nonlinear I-V behavior of the nano-device. How- 
ever, to be able to observe this effect, conducting channels 
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in the device need to be present, otherwise transmission 
saturation is reached at small voltages and a linear I-V 
characteristics is seen. Formally, it is the transmission 
coefficient T{E, Vb) that is of crucial importance when 
analyzing the effect of the eigenstates of the leads and 
the device, as well as the lead-device coupling. 

In the next section we determine the behavior of 
T{E, Vb) and we will address the following questions: Is 
the voltage independence of T{E, Vb) an adequate pic- 
ture? Which T{E, Vb) behavior would result in nonlin- 
ear I-V characteristics? Can the major characteristics of 
I{Vb) be qualitatively estimated from simple arguments 
or does one always need to perform an extensive ab-initio 
simulations? 



tively. Due to the conductance channel saturation ef- 
fect in the scattering region, the transmission coefficient 
has a very weak energy dependence (roughly constant). 
This behavior also emerges in the I-V curve which shows 
a very weak nonlinear behavior. We conclude that the 
roughly linear I-V curves of FigJ|-B and Figi-B,C are 
due to very different origins. In the latter case it is due 
to the channel saturation effect in the scattering region, 
whereas in the former it is due to a compensation be- 
tween the effects of increasing energy and bias voltage 
on T{E,Vb). 

V. DISCUSSIONS 



IV. BEHAVIOR OF TRANSMISSION 
COEFFICIENT T{E,Vb) 

For all the Au nano-contacts we have investigated the- 
oretically, T{E, Vb = 0) increases as a function of E (for 
E < 0.2eV). A typical behavior is shown in Figj^a by 
a dashed line for an ideal top site, pure and perfect Au 
device. From this curve, it is clear that there is trans- 
mission enhancement as a function of E which should 
result in a nonlinear I-V curve if T{E, Vb ^ 0) behaves 
in the same way. In the same graph we have also plotted 
T{E, Vb = 0.136^). One can observe that the general en- 
ergy dependence of E is still the same, namely increasing, 
but there is a global shift of the curve downward. There- 
fore, an increase in transmission coefficient as a function 
of E is compensated by its decrease due to increased bias 
voltage. The total effect on the current is to produce a 
linear I-V curve as shown in Fig.^B. We also note that 
this complete compensation between E and Vb is not uni- 
versal and can be different from one system to another. 
In fact, even for the same device it can behave differently 
at different energy ranges. This results in different fea- 
tures in the I-V curves such as plateaus, the "wiggles", 
etc. as we have discussed previously. 

A similar analysis is done on a device with a S impu- 
rity. The results are shown in Fig.^b. For this device 
it is clear that the effect of E on transmission at Vfc = 
is more pronounced. This is evidence that tunnelling is 
an important factor. In addition to this, we note that 
applying a bias voltage to the system causes a decrease 
in the transmission coefficient, but it is not a global de- 
crease. In particular, a bias voltage of Vb = 0.136V actu- 
ally increases transmission at small energies as shown in 
Fig.^b. Therefore, the combined effect of E and Vb does 
not cancel and it gives rise to the nonlinear I-V char- 
acteristics as seen in Fig.|^. According to this analysis, 
one can easily predict some aspects of the I-V curves just 
by studying the transmission coefficient T{E, Vb) at two 
different bias voltages. Obviously, this helps to predict 
I{Vb) characteristics with much less computational effort. 

In Figj^Cjd, we do a similar analysis on a device with 
disorder and one with disorder plus impurity, respec- 



We have already shown in the previous sections that 
transport at the molecular level is a complex phe- 
nomenon. To understand these systems, careful experi- 
mental work which separates electronic effects from struc- 
tural relaxations, as well as detailed calculations which 
include the effects of the molecule, leads and their cou- 
pling are required. In this section we discuss and com- 
pare our findings with previously published theoretical 
concepts and experimental results. 

Experimental work reported by Costa-Kramer et. 
oLEjEj have shown clear nonlinear I-V characteristics in 
Au nano-contacts starting at bias Vb = Q.IV, and its 
origin was attributed to strong e-e ir^ractions. To rule 
out the impurity effect, the authorsE^ have used scan- 
ning electron microscopy to analyze contacts of diameter 
^ SOOnm. They also used energy dispersive X-ray anal- 
ysis and determined a contamination concentration be- 
low detection sensitivity. Experimental cleanness checks 
performed for the large nano-contacts 300nm), how- 
ever, can not be extrapolated to junctions of few atoms 
in size due to the exquisite chemical sensitivitv-as demon- 
strated by our model. In other experimentsH with Au 
relays, it was observed that conductance quantization 
histogram survived for even larger bias, (1 x Go peak 
persists for Vb ~ 1.8^). These data show clearly that in 
such an experiment nonlinear I-V behavior cannot occur 
at low bias {Vb ~ O.iy). We believe that these junc- 
tions were formed between atomically clean gold con- 
tacts. These experiments were performed by forming 
and breaking the contacts very quickly (^ fJ'Sc), thus 
removing the mpurity atoms from the Au junction even 
if they existed lA. Recently, elegant experimental work by 
Hansen et. alzB found that contaminated Au nanostruc- 
ture show nonlinear I-V characteristics, whereas experi- 
ments done with clean tip-sample in UHV show perfect 
linearity for Vb < 0.7V. The nature of the contamination 
was not determined. 

A fundamental question that is to be addressed in this 
section, is transport through an impurity the most likely 
physical picture that can explain the observed nonlin- 
ear I-V characteristics. In the following we compare and 
discuss some of the concepts described in the literature 
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relevant to this issue. 

Free electron models have been used to describe the 
behavior of nanostructures under external biasEirEj. In 
these systems, depending on the potential profile across 
the device and the external bias voltage drop, various I- 
V characteristics can be extracted. In passing, we note 
that these calculations neglected the voltage dependent 
couphng T{E, Vb) (discussed in Fig.Q). This is obviously 
a major draw back for calculations at the molecular level, 
as shown by our analysis. However, a voltage drop at the 
contacts is still a reasonable approximation. In Figj^, 
we plot the Hartree potential across the four Au-atom 
chain. It is clearly seen that the potential drops mostly 
at the interfaces. However, assuming a uniform potential 
across the constriction is not adequate due to the atomic 
structure and the small variation of charge transfer as a 
function of bias voltage. 

A self consistent tight binding (TB) model has also 
been implemented to find the conduction dependence of 
each eigenchannel in Au nano-contact as a function of 
bias voltageE3. In these calculations, the TB parameters 
are calculated from a bulk system and charge neutrality 
of each atom of the nano-contact is enforced in order to 
carry out the self consistent calculations. Although it 
is not clear if the TB parameters determined from bulk 
structures are directly transferable to nano-contacts with 
atoms of low coordination number, especially when put 
under a bias potential, our ab initio results show that 
charge neutrality is a valid approximation for Au devices 
without impurities since charge transfer is quite small. In 
self consistent tight binding models, charge neutrality is 
accomplished by locally adjusting the chemical potential. 
For a zero bias calculation and nano-contacts with three 
atoms, Ref. ^ shows that a local potential ^ 3eV needs 
to be added to the central atom to achieve charge neu- 
trality. This seems to be a large value for Hamiltonian 
correction. Therefore, we suggest that for pure Au nanos- 
tructures, zero bias ab initio calculations should be done 
to extract the TB parameters (rather than from bulk). 
These can then be used more safely with charge neu- 
trality constraints to deduce TV characteristics of nano- 
contacts. Since we have seen only small effects of the bias 
voltage on the atomic charge transfer for a S doped Au 
structure in our ab initio calculations, we suggest that 
for this particular structure it is possible to deduce the 
charge distribution and the TB parameters from zero bias 
calculations. Using these TB interaction parameters and 
the zero bias charge at each atom as a constraint, the 
I-V characteristics of these structures can then be solved 
self-consistently. The results are less accurate compared 
to a full ab initio calculation, but they should give bet- 
ter results than using the conventional TB parameters 
derived from bulk systems. 

There are other important theoretical calculations go- 
ing beyond the single particle picture. Since the device 
constriction is narrow, e-e interaction can be strong and 
non-Fermi liquid behavior might have to be taken inta 
account. However, it was shown by Maslov and Stoneo 



that for these systems the resistance is due to contact 
and thus the results are independent from the Luttinger 
liquid behaKior in the constriction. Therefore, in these 
calculationsEj transmission is found to be identical to 
non-interacting particles. An important, single parti- 
cle transmission, approximation is incorporated in this 
modelEj. However, at non-zero temperature or/and bias 
voltage, a finite number of particles are injected into the 
nano structure resulting in backscattering effects. These 
lead to charge accumulation at the interface, creating an 
extra potential in addition to the original constriction po- 
tential. This additional potential is both Hartree (Vtr) 
and exchange-correlation (Vxc) in nature- and it is the 
reason for the so called resistance dipoldlJ. Ehie to this 
additional potential contribution, it was showiO that for 
ID systems the transmission coefficient is renormalized. 
It was further proposed that this charging effect can evem 
close a conducting channel that is 90% transmissiveo. 
This channel is gradually opened as Vb is increased for a 
complete transmission at Vb ^ 0.35V^, thereby inducing 
a nonlinear I-V curve. From a theory point of view, the 
picture of charging induced nonlinear I-V characteristics 
should overcome two further difficulties: that the renor- 
malized transmission depends on an interaction parame- 
ter a which can not yet be determined for atomic wires; 
and that previous calculations solved a ID case with no 
effect of Vb on the transmission T{E,Vb). In the rest of 
this section, we follow the interesting idea of the charg- 
ing effect and analyze it in more detail to .understand if 
this effect, which leads to channel closingE^EJ can give 
rise to nonlinear I-V curves of aimnic devices. To start, 
we follow the work of Yue et. alL3 by assuming a strong 
interaction_apd write the renormalized transmission co- 
efficient asHO, 



To{E/Do) 



Ro 



(4) 



To{E/Do)^'^ 

where, To,Rq are the transmission and reflection coeffi- 
cients of the non-interacting model such that Tq+Rq = 1, 
and a is a parameter to describe e-e interaction in the 
constriction. The parameter Dq is the energy range 
near Ep wliich contributes to renormalizing Tq and is 
determinedO by Dq = Tivp/W, where vp is the Fermi 
velocity of the system {lO^craJsc for Au) and W the 
width of the nano-constrictiorO, Vt^ 10 — 20A. There- 
fore, we find that for these devices Dq l.OeV. From 
Eq. ^, we compute the current I{Vb) at zero tempera- 
ture by integrating T^{E) from zero to eVb, assuming 
no bias dependence of T^{E). The conductance is then 
deduced by G = I{Vb)/Vb- We notCrthat due to the fi- 
nite length effect of the constrictiona, L ^ lOOA, the 
renormalization effect is cut off for bias voltages Vb < Vg, 
where Vs — 2'iThvp/L ^ O.IZJq- Within this approach 
and using the free parameters as specified, we have plot- 
ted in Fig.|^a the conductance of a single channel as a 
function of bias voltage Vb- Qualitatively, the results 
show that G increases with Vb due to the channel open- 
ing. However, quantitatively they do not give a complete 
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channel opening at the experimental value of Vb ^ 0.35V 
(as suggested in Ref .p^ , if the channel is less than 10% 
transmissive at V;, = 0. In fact, we found that a bias of 
2V is needed to overcome the charging potential barrier, 
thereby the nonlinearity in I-V curve can only set in at 
much larger voltages. We also note that since the extra 
charge due to backscattering is accumulated at the in- 
terface which has a larger cross section, we expect that 
DFT and the local density approximation to Vxc should 
work well. Our self-consistent calculations, with all the 
charge transfer and re-arrangements accounted for, have 
already partially included the backscattering effects. Our 
results show that it is possible to partially close a chan- 
nel ^ 60%-but not completely. Our experimental data 
which shows zero bias conductance of ~ 2.2Go, indicate 
that at least one channel is 20% transmissive and that 
channels are only partially closed, consistent with our the- 
oretical work. Therefore, we can conclude at this point 
that the physical picture of electron interactions to com- 
pletely close and open a channel as a function of Vb, while 
interesting, cannot explain our experimental data. 

To fiaxther address the importance of backscattering 
effectsEil, especially for our devices with realistic atomic 
leads, we have investigated the dependence on the inter- 
action parameter a in Eq.(Q). Tuning this parameter can 
dramatically change the renormalized transmission coef- 
ficient as a function of Vb as shown in Fig.^b. There- 
fore we have a wide range of this parameter for fitting 
our experimental data. However, if we change the Au 
nano-contact by only a single sulphur atom, we should 
not expect a very different e-e interaction parameter a, 
thereby this would predict an I-V curve not too different 
from the one without the sulphur. This suggests that 
strong interaction alone would not predict strong sensi- 
tivity of the I-V curve dependence on small amou pt. . o f 
impurities, in contradiction to experimental resultsP^IM. 
Transmission renormalization is important in ID scatter- 
ing problems where there is no charge in the leads. For 
realistic atomic leads, the backscattered charge is a small 
fraction of the original one, and its effect should be well 
screened. With all these considerations, we believe it is 
unlikely that the dynamically established charging effect 
alone is large enough to cause the observed I-V nonlin- 
earity for Au nano- contacts. 

VI. CONCLUSION 

In this paper, we discuss the electronic transport prop- 
erties of Au nano-contacts from both experimental and 
theoretical studies. Our experimental data analysis en- 
ables us to separate electronic effects from structural re- 
laxations, allowing a better comparison to theoretical 
modeling of these nano-devices. Our theoretical work 
shows that transport properties at the molecular scale 
need to be analyzed at the systems level: leads, molecule 
and their interactions have to be studied simultaneously. 



The self-consistently determined transmission coefficient 
T{E,Vb) is shown to vary as a function of E and Vb- 
This gives rise to differential conductance fluctuations of 
the order of 1% at temperature 300K taking into ac- 
count the experimental averaging process. These fluc- 
tuations are attributed mainly to the effects of the lead 
band structure. Most striking, however, is the fact that 
a single impurity atom at the contact region can alter 
I-V curves qualitatively in these devices: pure and per- 
fect Au nano-contacts do not give the observed nonlin- 
earity, while a sulfur doped device does. Importantly, 
we have shown that the measured nonlinear I-V charac- 
teristics of Au nanowires can be quantitatively modeled 
by impurity effects which create a tunnelling barrier at 
the nanostructure junction. This effect, among others, 
points to the vital importance for understanding contacts 
in nano-electronic devices, and, perhaps, to exploit it for 
the bencflt of device operation. 
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FIG. 1. (a)-Schematic of the experimental set-up with 
a mechanical controllable break junction (MCBJ) shown at 
the centre with An wires attached. (b)-DC conductance I /V 
(lower) and differential conductance SI/5V (upper) plotted 
vs V; Inset is the I-V curve. The dashed line, with slope 
~ 2.2Go, corresponds to linear behavior and helps to rec- 
ognize nonlinearity. These experimental results are measured 
over a 5 second voltage sweep from positive (dark) to negative 
(grey) bias. 

FIG. 2. (a)-The Shape function (5',Eqj|) vs V calculated 
from the data shown in Figj^b. (b)-Normalized functional 
form of the voltage fo{V) with linear dashed line shown to 
help view the onset of non linearity, (c)- Unbiased Conduc- 
tance g{X) vs V. (d)-Normalised conductance, 5f/5V vs V; 
with the inset corresponding to a magnified section of the 
same data to show more clearly fine details. In all graphs, 
dark (grey) lines correspond to positive (negative) bias. 



FIG. 6. I-V characteristics of Au contacts with 2 S impurity 
(A), disordered interface (B) and disordered interface with 1 
S impurity (C). The insets illustrate the corresponding struc- 
tures, with S represented as a dark circle; Au atoms in the 
lead and the chain are shown with grey and light grey circles, 
respectively. (A), (B) and (C) have zero-bias conductance of 
O.SSGo, 0.915Go and 0.833Go, respectively. 

FIG. 7. Scaled transmission coefficient T{E) at Vi, = 
(dashed line) and 14 — 0.136V (solid line) for ideal contact 
structure at the top site (a), IS impurity structure (b), disor- 
dered interface (c) and disordered interface with 1 S impurity 
(d). The scaling factor is T{E = 0, 14 = 0). 

FIG. 8. (a) Average Hartree potential (Vh) across the cen- 
tral cross section (~ 9A^) of the ideal contact at the top site; 
dotted line at zero bias, and dashed line at Vb = 0.5V; and 
(b) their difference showing that most of the potential drops 
symmetrically at the interfaces. Black circles correspond to 
atom positions along the constriction for guidance 



FIG. 3. TV characteristics of Au contacts with Hollow site 
registry (A) and Top site registry (B); with zero- bias con- 
ductance 0.94Go, O.SGo for (A) and (B), respectively. The 
structure of the device is illustrated in (a) where LL, D and 
RL correspond to Left Lead, effective Device and Right Lead, 
respectively. The Hollow site is shown in (c) and Top site (d) 
where chain atoms are illustrated by dark circles. Inset (e) 
shows the band structure of the lead along the transport di- 
rection, the conducting band is shown by solid line. Inset (b) 
shows differential conductance fluctuations of the Hollow site 
as a grey dotted line, dark dotted line and dark continuous 
line for T = OK,T = 300K and for T = 300K with 4 mV™^ 
modulation voltage taken into account; respectively. 



FIG. 9. Conductance calculated from renormalized trans- 
mission due to backscattering as a function of bias voltage; (a) 
for To = 0.99, 0.9, 0.8 corresponding to solid, dashed and dot- 
ted line, respectively with interaction parameter q = 1.0 and 
(b) for To — 0.9 and interaction parameter a — 0.1,0.5, 1.0 
corresponding to solid, dashed and dotted line, respectively. 



FIG. 4. Contour plots of surface charge density at par- 
ticular bias voltage Vb{V) and Energy E{eV) indicated on 
each figure. We use the same scale for all graphs to compare 
the conductance probability for each configuration. Dark cir- 
cles correspond to atom positions along the constriction for 
guidance. 



FIG. 5. TV characteristics of Au contact doped with S im- 
purity (Solid line). The dashed line has a slope ~ 0.67Go 
and is shown to help view the onset of non linearity; Cir- 
cles correspond to the experimental results obtained from 
fo{V) (Fig^-b) and multiplied by the zero bias conductance, 
0.67Go. The inset (a), shows the atoms registry, where dark 
circles show the atomic chain, grey ones are Au atoms in the 
leads' surface and light grey is the S atom. 
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